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Abstract 


Voyager images of Triton indicate considerable spatial variability in the concentration of at 
least two different scattering components in the atmosphere. Data from high phase angle limb 
scans have been fit to IVlie scattering models to derive mean particle sizes, number densities, and 
vertical extent for both types of scattering material at ten different locations between 15 S and 
70°S. These fits reveal a thin haze at latitudes equatorward of 25-30°S. The imaging data can be 
fit reasonably well by both conservatively scattering and absorbing hazes with particle sizes near 
0.18 \l m and optical depths of order 0.001-0.01. Rayleigh scattering haze fits the imaging data 
somewhat less well, and can be totally ruled out by combining the imaging and UVS measure- 
ments. At high southern latitudes, Triton displays clouds below an altitude of ~8 km, as well as 
the haze at higher altitudes. The clouds have particle sizes which may range from 0.7-2.0 fim, or 
may be near 0.25 fjxti. The atmospheric optical depth poleward of 30 S must be generally greater 
than 0.1, but need not be more than 0.3. Horizontal inhomogeneities are quite noticeable, espe- 
cially at longitudes east of (i.e. higher than) 180°. 
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1. Introduction 

Until the Voyager encounter with Neptune very little was known about its largest satellite Triton. 
Even estimates of Triton’s radius were uncertain. It was known from measurements of the 
infrared spectrum (Cruikshank et al. 1984, Cruikshank etal. 1989) that condensed nitrogen was 
present in some form on the surface, but the surface pressure and temperature and the atmospheric 
optical depth were unknown. From radio occultation data taken during the Voyager encounter 
(Tyler et al. 1989) it was determined that Triton possessed an atmosphere which was mostly nitro- 
gen, with a surface pressure of about 1 2 /ibar and a surface temperature of 37K. Despite this ex- 
ceedingly low pressure, images of Triton’s limb taken near closest approach clearly show some 
kind of scattering layer above the surface (Smith et al. 1989). A detailed analysis of these images 
(Pollack et al. 1990) indicated that a thin haze, presumeably of photochemical origin, is present 
nearly everywhere on Triton. In addition, a sizable fraction of the surface poleward of 30 S is 
covered by discrete clouds which are confined to the lower 8 km of the atmosphere and may have 
optical depths of order 0.1. The confinement of these clouds to the lower part of Triton’s atmo- 
sphere and high southern latitudes suggests that they are composed of nitrogen ice. 

In this paper, we will present the results of a detailed analysis of the Voyager images of 

Triton taken at the highest solar phase angles. Although the spatial resolution of these images is 

lower than that of the images studied in Pollack et al. (1990), observations at high phase angles are 

more sensitive to properties such as the particle sizes and vertical extinction profiles of the hazes 

and clouds than the highest resolution images were. In the next section, we outline our data anal- 
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ysis procedure. In sections >2?*and.3^we discuss the properties of atmospheric models derived for 
low and high southern latitudes, respectively. SectionyTcon tains a more detailed discussion of lo- 
calized atmospheric features, and the final section considers some implications of our findings for 
models of seasonal changes on Triton. 
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2. Data set 


The data set consists of fifteen Voyager wide angle images of Triton taken at high solar phase 
angles (140°-160°). Three of these images were taken with the violet filter, five with the blue fil- 
ter, five with the green filter, and two with the orange filter. The linear resolution in these images 
ranged from about 10 kilometers per line pair (km/lp) at the beginning of the sequence (highest 
phase angle) to about 25 km/lp at the end of the sequence (lowest phase angle). Details are given 
in Table 1. 

Initial processing of these images was the same as that outlined in Pollack et al. (1987), and 
included removal of blemishes and reseaux, dark current subtraction, and conversion (via a non- 
linear response function specific to each pixel) from DN value (0-255) to specific intensity (I/F) at 
each pixel. The precise pointing of each image was initially determined by fitting Triton’s ob- 
served limb in the image to a circle of radius 1350 km, using the NAV program in the VICAR 
image handling package. Later the pointing was improved by comparing specific intensity profiles 
calculated from the best fitting atmospheric models to the observed I/F scans and determining 
whether all the limb scans for a given image could be better fit by shifting the observed limb. 

Specific intensity was measured as a function of line-of-sight radius at intervals of approxi- 
mately one pixel along selected radial lines crossing the limb of each image. The scan lines for all 
the images are shown projected onto the surface of Triton in Figure 1. As can be seen in this 
figure, the scan lines were chosen to correspond to ten distinct areas along the limb. The numerical 
analysis described in the following sections was carried out by using all the scans comprising a 
particular region to form a single data set covering a range of solar phase angles and wavelengths. 

Before fitting the limb scans to the atmospheric models a final background subtraction was 
carried out by fitting a quadratic function to the portions of each limb scan well above the limb and 
well inside the terminator. This is necessary at high phase angles because some sunlight is scat- 
tered from the sunshade into the camera, and this is not included in the standard dark current sub- 
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traction. Also, the dark current is proportional to the temperature, and it is generally not possible 
to find a dark current frame taken at the same temperature as a given image. 

3. Low latitudes (haze only) 

The specific intensity as a function of radius was calculated for an atmosphere consisting of con- 
centric homogeneous spherical shells using methods described in Rages and Pollack (1983) and 
Pollack et al. (1987). The spatial resolution of the high phase angle Triton images was not high 
enough to obtain meaningful vertical extinction profiles by directly inverting the radial intensity 
profiles. Instead, the radial intensity profiles were fit to simple exponential models of the extinc- 
tion as a function of altitude. 

At latitudes near the equator the atmospheric model contained a single component (haze), 
with a number density varying exponentially with altitude. The single scattering properties of the 
haze were described by Mie theory, with a complex refractive index and particle size that did not 
vary with altitude. The real part of the haze refractive index was assumed to be 1.285, consistent 

with the value for methane ice. The imaginary part of the refractive index was assumed to be a 
linear function of wavelength, n f =a+bX . The particle sizes were distributed according to a fairly 

narrow (width parameter 0.05) Hansen-Hovenier size distribution with a cross section weighted 
mean radius r h . 

It was assumed that there were no horizontal variations in the atmospheric properties over the 
entire area covered by the limb scans in any one region. This is a reasonably good assumption at 
low latitudes, but some of the higher latitude regions cover more than 30° of latitude, and in some 
cases (discussed in more detail in section >*5 this assumption is clearly violated. 

Triton’s surface photometric properties were described by a single Hapke photometric func- 
tion for each filter, with parameters given by Hillier et al, (1991) for the violet, blue, and green fil- 
ters. Orange filter parameters were estimated by extrapolating from blue and green. 
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Specific intensities were calculated at 30-50 points along each radial scan at ten wavelengths 
between 0.3 /im and 0.65 /im, then convolved over both the various filter response functions 
(Danielson et al. 1981) and the point spread function (S.A. Collins, private communication) of the 
Voyager camera. 

The haze column number density n ^ , scale height H ^ , mean particle radius r ^ , and wave- 
length-dependent imaginary part of the refractive index were treated as free parameters in a least 
squares fitting procedure. 

Results are shown in Figure 2 for a region near 15°S, 275°E (region 10). Three separate 
scattering models were found which gave acceptable fits to the imaging data: Rayleigh scattering, 
conservatively scattering particles with a mean radius of 0.19 fjm, and 0.18 /im particles with an 
imaginary index which is moderately large in the violet and decreases toward longer wavelengths. 
The details of these three models can be found in Table 2, together with divided by the number 
of degrees of freedom in each data set. The (radial) scattering optical depths for these models gen- 
erally lie between 0.001 and 0.01 in all four filters, as shown in Table 3. 

The Rayleigh scattering model fits the imaging data somewhat less well than the other two, 
and can be eliminated entirely by a comparison with data from the Voyager Ultraviolet Spectrome- 
ter (see Krasnopolsky et al. 1992). The Rayleigh scattering model predicts ultraviolet (slant) opti- 
cal depths to the surface of around 10, which is an order of magnitude larger than those seen by the 
UVS. The other two models predict UV optical depths of about 0.15 if the imaginary index in the 
UV is given by extrapolating the linear function used in the visible. In fact, photochemically pro- 
duced hydrocarbons could well have an absorption edge between 0.3 /im and 0.165 /im (the long 
wavelength limit of the UVS data) which could raise the UV optical depths to the 0.5- 1.0 range 
observed. 

The data from regions 1 and 9 were also fit to single component scattering models, and the 
details of these models are also given in Table 2. In both cases the derived imaginary index is 
very steeply sloping (in one case rising and in the other falling toward the red). 
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4. High latitudes (haze + cloud) 

In light of the findings of Pollack et al. (1990) that both the haze and optically thicker discrete 
clouds can be present south of about 30°S, at these latitudes a two component model was fit to the 
data. This model contained the 0.19 jim conservatively scattering haze above an altitude of 8 km, 
with the column number density n ^ and scale height left as free parameters. Below 8 km there 

was another scattering component (cloud), with an independent exponential altitude profile defined 
by n Q and H c , particle radius r. , and wavelength dependent imaginary part of the refractive index. 

The real part of the cloud refractive index was taken to be 1.25, consistent with measured values 

for nitrogen ice (Thompson et al. 1973). 

This model was applied to all the regions shown in Figure 1 except region 10. It soon 
became clear that large cloud optical depths (-1) were favored by the fitting algorithm. Optical 
depths this large are excluded by the images taken at lower phase angles and higher spatial resolu- 
tions, which show little or no indication of surface obscuration. Therefore an additional constraint 
was included in the fitting algorithm, specifying an upper bound Tj, on the total atmospheric optical 


depth. The function to be minimized was then 
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x 2 = 
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n =1 
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where the first summation is the usual formulation of % 2 as the weighted sum over all N data 
points of the squares of the differences between the observed and calculated values of the depen- 
dent variable (in this case specific intensity), and 
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where <r T controls the “hardness” of the upper bound. Initially was set to 0.3 and <7 r was set to 
0.05. 

Results for a typical high latitude region (region 3, near 50 S, 124 E) are shown in Figure 
3. There are two distinct scattering models which fit the data for this region almost equally well. 
One contains cloud particles with a mean radius of 1.5 /rm (solid line), while for the other this 
quantity is 0.22 jum (chain dotted line). Figure 4 illustrates how these two very different particle 
sizes could both fit the observations. It shows the albedo-weighted single scattering phase func- 
tion for both models and all four filters. Although the larger particles have a much larger diffrac- 
tion peak than the smaller particles, it is confined to scattering angles less than 20° (phase angles 
>160°), outside the limi ts of the Voyager observations. Between 20° and 40° scattering angle, the 
range which is covered by the observations, the two sets of phase functions are very similar. 

The effect of lowering the atmospheric optical depth on the fit to the data was found by 
changing Tj and <7^ to 0.08 and 0.02 respectively, which limited the optical depth to around 0.1. 

The result is shown as the dashed line in Figure 3. Forcing the optical depth to be this low in this 
region increases % by a factor of 5. 

The cloud scale height is more than 16 km, making the total thickness of the cloud layer less 
than half a scale height. It was not possible to determine whether the cloud was uniformly mixed 
with the gas, which has a scale height of around 18 km (Tyler et al. 1989), or whether the cloud 
mixing ratio increased with altitude as might be the case if condensation were occuring at the top of 
the cloud layer. 

The best fitting parameter values for each two-component model studied for all of the regions 
1-9 are given in Table 4. The total (vertical) optical depths for all the models (both one- and two- 
component) are given in Table 3. Most of the high latitude regions were similar to region 3 in 
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having two distinct particle size regimes that fit the data almost equally well. At the highest lati- 
tudes (regions 5, 6, and 7) r m was near 0.7 jim or 0.3 /im. Region 9 was greatly affected by hori- 
zontal inhomgeneities in the atmospheric scatterers (see next section). Haze column number densi- 
ties at high latitudes are about the same as they are at low latitudes (although at high latitudes the 
haze optical depth is given for 8 km altitude, while at low latitudes it is given at the surface). The 
two regions for which both one- and two-component models were tried (regions 1 and 9) the two- 
component models give significantly better fits to the observations. 

5. Horizontal inhomogeneities 

Many regions (especially those east of 180° longitude) showed features in the radial intensity pro- 
files that can only be due to local variation in the atmospheric scattering properties. Figure 5 
shows the complete data set from region 6, near 70°S, 175°E. In this region, most of the limb 
scans show a pronounced hump on the portion of the curve inside Triton’s limb (radius less than 
1350 km), represented by open circles in Figure 5. The surface is not bright enough in the scatter- 
ing geometries represented by these scans for these humps to be due to variations in the surface 
properties, and this behavior cannot be produced by any extinction profile composed of homoge- 
neous spherical shells. The most likely cause of the humps in this data set is a local increase in the 
cloud optical depth over part of the area covered by region 6. This is a strong indication that the 
“background” cloud optical depth in this region must be considerably less than unity. 

The data points shown as open circles in Figure 5, and similar data points in the other re- 
gions, were not included in the data set used to define the values of the free parameters given in 
Table 4, because the two-component atmospheric model itself is inadequate for reproducing them. 
The location of all such data points are shown in Figure 6. 

The most pronounced local variations in Triton’s atmospheric scatterers seem to be concen- 
trated in a diagonal band between 55° and 15° south latitude, and 180° and 270° east longitude. 
A comparison of Figure 6 with Figure 1 shows that this apparent concentration is not merely a se- 
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lection effect. The pan of Triton covered by the high phase angle images is (more or less in- 
evitably) on the opposite side of the satellite from the region covered at the highest spatial resolu- 
tion at lower phase angles. Therefore the characteristics of the surface in the band of local atmo- 
spheric variability are known only vaguely. According to McEwen (1990), this band lies mostly 
over the bright south polar cap, with its northern end extending into the brighter polar cap fringe. 
The band of pronounced atmospheric variability lies along the evening terminator, so it could be an 
indicator of some kind of diurnal effect. 

The plumes visible on the highest resolution images of Triton, marked by x’s in Figure 6, 
are also in the opposite hemisphere from the areas of local variability seen at high phase angles. 
One area which Pollack et al. (1990) noted as unusual because it lacked the otherwise ubiquitous 
haze is somewhat closer to the variable regions (dashed ellipse in Figure 6). 

6. Discussion 

Sagan and Chyba (1990) give an expression for the sedimentation velocity of particles in Triton’s 
atmosphere as a function of particle size. Assuming a surface pressure of 12 ,ubar, and an 
“equivalent isothermal temperature” of 48K (from Tyler et al. 1989) to characterize the variation of 
pressure with altitude, this expression gives fall rimes to the surface from an altitude of 8 km of 
order 1 day for 1.5 }Jm particles and 1 week for 0.23 jjm particles. 

Using the coagulation kemal given by Turco et al. (1982), simplified by assuming mon- 
odesperse particle sizes, coagulation times are found to be of order 1 month for small cloud parti- 
cles and 1 year for large ones. Condensation times are critically dependent on the difference be- 
tween ambient pressure and the saturation vapor pressure, neither of which is known well enough 
to yield useful limits. 

The fact that the coagulation times are substantially longer than the fallout times may consti- 
tute further evidence that condensation is taking place, but it is not possible to discriminate between 
the two possible cloud particle size regimes by comparing growth and fallout times. 
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Nitrogen is the major constituent of Triton’s atmosphere, and its ambient pressure is essen- 
tially the same everywhere on Triton. If the clouds are composed of condensed N 2 , the fact that 

they appear preferentially toward the south pole would indicate systematically lower atmospheric 
temperatures there, despite the fact that solar insolation is currently highest at high southern lati- 
tudes. Spencer and Moore (1992) have postulated that the high southern latitudes of Triton may be 
mostly frost free, with the high surface albedo produced by a non-volatile substrate such as f^O or 

C0 2 ice, and at least some of their frost-free models have surface temperatures which decrease 

monotonically between the equator and 60-70°S. To the extent that surface temperatures on Triton 
can be correlated with atmospheric temperatures, this would be consistent with increased N 2 con- 
densation at higher southern latitudes. 
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Table 1: Description of Triton Images 


FDS count 

Camera 
and Filter 

Phase angle 
(°) 

Resolution 

(km/lp) 

ll39~8.h 

WA-B 

159.2 

9.5 

11398.17 

WA-G 

157.8 

10.1 

11398.23 

WA-V 

156.4 

10.6 

11398.29 

WA-B 

155.0 

11.2 

11398.35 

WA-G 

153.5 

11.7 

11398.56 

WA-V 

150.2 

12.7 

11399.07 

WA-G 

148.0 

14.8 

11399.18 

WA-B 


15.9 

11399.29 

WA-0 

146.1 

17.0 

11399.51 

WA-B 

144.5 

19.1 

11400.02 

WA-G 


20.2 

11400.13 

WA-V 


21.4 

1 1400.24 

WA-B 

143.8 

22.4 

11400.48 

WA-0 

141.3 

24.8 

11400.54 

WA-G 

141.3 

25.4 
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Table 2: Optimum Values of Free Parameters for 1-Component Models 


Region 

model 

r h (pm) 

n h 

(10 8 cm* 2 ) 

H h 

(km) 

a 

b 

(pm 1 ) 

x 2 

/d.o.f. 

1 

a) 

0.1 6210.007 

0.22±0.05 

5.4±0.3 

-0. 5+0.4 

1.6±0.8 

0.80 

9 

a) 

0.0962±0.0012 

3.4±0.5 

5.1±0.3 

1.1±0.3 

-1.6±0.4 

1.28 

10 

a) 

0.182±0.017 

0.037±0.018 

10.8±0.7 

0.12±0. 15 

-0.210.4 

0.79 


b) 

0.188+0.014 

0.029+0.010 

11.0+0.6 

0.0 

0.0 

0.80 


c) 

0.025 

3720±140 

10.8±0.6 

0.0 

0.0 

0.92 
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Table 3: Model Optical Depths 

Vertical optical depth to surface 


Region model Violet Blue Green Orange 


1 

a) 

0.027 

0.028 

0.031 

0.032 


b) 

0.127 

0.126 

0.125 

0.125 


c) 

0.039 

0.040 

0.042 

0.041 

2 

a) 

0.309 

0.307 

0.304 

0.308 


b) 

0.364 

0.319 

0.237 

0.208 


c) 

0.111 

0.100 

0.082 

0.076 

3 

a) 

0.364 

0.360 

0.357 

0.361 


b) 

0.409 

0.353 

0.254 

0.221 


c) 

0.117 

0.104 

0.084 

0.078 

4 

a) 

0.369 

0.364 

0.361 

0.366 


b) 

0.411 

0.355 

0.255 

0.221 


c) 

0.121 

0.104 

0.077 

0.068 

5 

a) 

0.351 

0.358 

0.352 

0.345 


b) 

0.374 

0.346 

0.294 

0.276 


c) 

0.112 

0.102 

0.086 

0.080 

6 

a) 

0.362 

0.355 

0.328 

0.317 


b) 

0.360 

0.345 

0.328 

0.324 


c) 

0.105 

0.100 

0.100 

0.101 

7 

a) 

0.358 

0.358 

0.342 

0.333 


b) 

0.377 

0.356 

0.314 

0.298 


c) 

0.109 

0.103 

0.095 

0.093 

8 

a) 

0.343 

0.344 

0.343 

0.342 


b) 

0.386 

0.323 

0.188 

0.136 


c) 

0.097 

0.096 

0.093 

0.091 

9 

a) 

0.143 

0.117 

0.065 

0.046 


b) 

0.151 

0.151 

0.151 

0.152 


c) 

0.331 

0.266 

0.164 

0.134 


d) 

0.089 

0.089 

0.087 

0.084 

10 

a) 

0.0047 

0.0039 

0.0024 

0.0020 


b) 

0.0036 

0.0030 

0.0020 

0.0017 


c) 

0.0143 

0.0100 

0.0050 

0.0039 
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Table 4: Optimum Values of Free Parameters for 2-Component Models 
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0.174±0.006 50±4 16±26 4.310.5 10.6±1.2 0.1±0.4 0.810.9 2.67 

2.01910.004 0.5210.04 16124 2.3810.15 11.711.0 -0.08510.013 0.2510 03 0 79 

0.0810.02 360015000 16121 2.6710.14 11.410.8 0.1810.17 -0.210.2 0.57 

0.10010.003 203114 16124 2.5410.16 11.210.9 -0.010.3 1.010.6 0.72 



Figure captions 

1. This figure shows the locations of the data sets fit to the spherical shell atmospheric models. 
The location of Triton’s limb in each of the green filter images is shown by the heavy lines, 
with the corresponding FDS numbers given in the legend. Each sheaf of thinner curved lines 
perpendicular to the limbs represents the radial limb scans taken from each of the 15 images, 
for each of ten labelled regions. These labels correspond to the region numbers in Tables 2, 
3, and 4. Finally, the location of the surface terminator is shown (thin solid line), as well as 
the location where the sun drops below the horizon for an altitude of 10 km above the surface 
(thin dashed line). Note that for some regions the limb scans cover more than 30° in latitude, 
which makes it rather unlikely that the atmosphere is homogeneous throughout the region. 

2. This is the complete data set for region 10, shown together with three possible model atmo- 
spheres: conservative scatterers with a mean particle radius of 0.19 jim (solid line), highly 
absorbing particles with a mean radius of 0.18 /im (chain-dotted line), and Rayleigh scatter- 
ers (dashed line). 

3. This is the complete data set for region 3. Three model atmospheres are shown: a) =0.3 opti- 
cal depths of 1.5 }im particles (solid line), b) =0.3 optical depths of 0.22 jj.m particles (chain 
dotted line), and c) =0. 1 optical depths of 0.20 fim. particles (dashed line). Model c) fits the 
data much less well than the other two. 

4. Albedo weighted single scattering phase function between 0° and 50° scattering angle, for 
models a) (solid line) and b) (dashed line) from Figure 3. 

5. This is the complete data set for region 6. This is an example of an “ill behaved’ 1 region, 
since many of the limb scans contain humps at radii less that 1350 km, which cannot possi- 
bly be fit by any horizontally homogeneous atmosphere, even if the condition of an exponen- 
tial extinction profile is relaxed. The points shown as open circles were not included in the 
least-squares fitting procedure. 
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6. The location of all the points in all the radial scans for which I/F was clearly affected by local 
inhomogeneities (most likely local increases in the optical depth) in the cloud properties. The 
size of each dot crudely corresponds to the amount by which I/F deviated from its 
“background” value. The size of two plumes seen in the highest resolution images of the 
surface are marked by x’s, and the location of a haze free area noted by Pollack et al. (1990) 

is shown as a dashed ellipsoid. 


19 



Location of limb scans 



100 150 200 250 300 360 

East Longitude 




0.0 0.1 0.2 0.0 




1400 1300 1200 

Equatorial Radius (km) 


o 

b 


1400 1300 1200 

Equatorial Radius (km) 



- /L 7‘ 























90 





















Triton Phase Function — 50°S Latitude 

1.48 /xm spherical haze 

0.22^tm spherical haze 



CNI 0 10 20 30 40 50 (NO 10 20 30 40 50 


Scattering Angle (°) Scattering Angle (°) 



(Nl 0 10 20 30 40 50 (NO 10 20 30 40 50 


Scattering Angle (°) Scattering Angie (°) 







1.0 




















East Longitude 




